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ABSTRACT 



A treatment is given of the orbit dynamics for linear unstable motion 
that allows for the zeros in the beta function and makes no assumptions 
about the realness of the betatron and phase functions. The phase shift per 
turn is shown to be related to the beta function and the number of zeros 
the beta function goes through per turn. The solutions of the equations of 
motion are found in terms of the beta function. 

1. Introduction 

In the case of linear unstable motion, the beta function can be zero at some points in the 
lattice. Because of the zeros in the beta function, and other assumptions often made about 
the realness of the beta function and phase function, the usual treatment given for stable 
motion does not carry over to the case of unstable motion. A treatment is given below, 
that allows for the zeros in the beta functions and does not make assumptions about the 
realness of the betatron and phase functions. 

It will be shown that the solutions of the equations of motion can be written in the 
form 



N z is the number of times /3(s) goes through zero between so and s. P indicates the 
principle value of the integral. The solutions of the equations of motion can also be 




(1-1) 
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written as 

x = exp [±ns/L\ f (s) (1-2) 

where f(s) is periodic and L is the length of one turn. It will be shown that for unstable 
motion 

H = 2tt (g + iq/2) 

q = \ Nz (1-3) 
" L ds 



p ri 



where N z is the number of zeros the beta function goes through in one turn. P indicates 
the principle value of the integral. 

Often, the case of unstable linear motion is found when a gradient perturbation is 
applied to a lattice whose unperturbed z/-value is close to q/2, q being some integer. In 
this case, perturbation theory will show [1] that the solutions have the form given by Eq. 
(1-1) where q/2 is the half integer close to the unperturbed to the z/-value. In the general 
case, where the unstable motion cannot be viewed as due to a perturbing gradient then 
the value of q is given by \N Z where N z is the number of zeros in the beta function in one 
turn. 

It will also be shown that near a zero of the beta function at s = si, ip will become 
infinite and the dominant term is ip is given by 

V>~±-log(s-si) (1-4) 



2. The Definition of the Beta Function 

The linear parameters can be defined in terms of the elements of the one period transfer 
matrix. The 2x2 transfer matrix, M, is defined by 

x(s) = M (s, so) x (so) 

-U) 

The one period transfer matrix is defined by 

M(s) = M(s + L,s) (2-2) 
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where the lattice is assumed to be periodic with the period L. The matrix M is assumed 
to be symplectic 

MM = 1 



M =SM S ( 2 3) 

-l o j V° v 

S is the transpose of S. Also \M\ = 1 where \M\ is the determinant of M. One can show 
that M(s) and M(so) are related by 

M(s)=M (s, s ) M (s ) M (s , s) (2-4) 

It follows from Eq. (2-4) that M\\ + M22, the trace of M, is independent of s. For 
unstable motion it is assumed that \M\\ + M22I > 2. This may be shown to lead to 
unstable exponentially growing motion. 

One can now introduce the constant parameter fx defined by 

cosh fi = - (Mh + M 22 ) (2-5) 

If Mn + M22 is positive, then fx will be real. However if M\\ + M22 is negative then fi has 
to have the imaginary part iqn where q is an odd integer. In general, one can write 

fj = IU R + i qn 

Mn + M 22 



cosh VR=7^ 



(2-6) 



where q is an even integer if Mn + M22 is positive, and q is an odd integer when Mn + M22 
is negative. It will be seen below that fiR is the exponential growth per period; that is, 
the growth per period for the unstable solution is exp( / u^). It will also be seen below that 
q is related to the number of zeros in the beta function, f3(s), in a period, which is 2q. 

fx is related to the eigenvalues of M, Ai and A2, where Ai + A2 = Mn + M22 and 
AiA 2 = 1 from |M - A/| = 0. It follows from Eq. (2-5) that 

Ai = exp (fj) 

(2-7) 

A 2 = exp (-fi) 

One can define the linear parameters, (3, a, 7, using the elements of the one period 
transfer matrix. If one uses the form of the transfer matrix often used [2] for stable motion 
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the linear parameters will be imaginary for unstable motion. To make the linear parameters 
real, they will be defined here in terms of the one period transfer matrix as 



M 



cosh [i + a sinh [i (3 sinh \i 

7 sinh fx cosh [i — a sinh /i 



(2-8) 



/3 7 = l - a 2 

f3, a, 7 are then given in terms of M y - as 

f3= (-l) q M 12 / sinh ja R 

a = (-l) q (m u - M 22 ) /2sinh^ (2-9) 
7 = (1 - a 2 ) /(3 

Eq. (2-6) does not specify the sign of One can define the sign of to be always 
positive. Then f3, a, 7 can then be computed from the My using Eq. (2-9). It will be 
seen later that the sign of (3(s) can change within a period, and (3(s) can be zero at certain 
values of s for unstable motion. 

Having defined the linear parameters, one can now find the relationships among them, 
their connection with the growth rate, the emittance and the solutions of the equations of 
motion. The treatment usually given for stable motion does not carry over for unstable 
motion, because if often assumes that (3 and ip, as defined for stable motion, are real and 
that (3 is never zero. One needs a treatment which does not make assumptions about the 
realness of /?, ip, and allows f3 to go through zero. This is given below. 

2.1. Differential Equations for /?, a, 7 

It is assumed that the linearized equations of motion can be written as 

dx . . 
— = A n x + A\2Vx 
as 

dp x , a (2-10) 

as 

An + A22 = 

In the large accelerator approximation, An = A22 = and A12 = 1. We note that 

-^-M(s,s ) = A M(s,s ) 

ds (2 _ n) 
— M(s ,s) = -M(s ,s)A 

CIS 
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The last equation follows from M(s, so)M(so, s) = I. Then using Eq. (2-4) 

M(s) = M(s,s )M(s )M(s ,s) 

one finds 



dM 
ds 



= A M - M A. 



(2-12) 



(2-13) 



A is the 2x2 matrix whose elements are the Aij of Eq. (2-10). Replacing M, using Eq. 
(2-8), in Eq. (2-13) gives the result 



^ = 2A n P - 2A 12 a 
ds 

da 



ds 



-A 21 p + A l2l 



(2-14) 



= 2A 21 a-2A lxl 

ds 

The first equation in Eqs. (2-14) gives the connection between a and f3 

1 ( ld P AO 

a= — ---f + A n p 
A 12 V 2 ds 



(2-15) 



2.2. Differential Equation for (3 

In this section, the differential equation for (3 will be obtained without making any assump- 
tions about the form of the solutions of the equations of motion. For the sake of simplicity, 
the derivation will be given for the large accelerator case which assumes An = A 22 = 
and Ai 2 = 1. 

Introducing 6, where f3 = 6 2 , Eqs. (2-14) can be written as 

db 

b— = —a 

ds 

da j 9 

dJ = A6+7 (2-16) 

? = -me 

ds 

K = -A 21 

The first two equations in (2-16) then give 



4- ( bf) = -Kb" 
ds \ ds 



= -Kb 2 



7 



b' 1 



(2-17) 



ds J 
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one then gets 

J? + Kb+ h\ (2-18) 

b = (3? 

Eq. (2-8) differs frm the usual result for stable motion only in the sign of the 1/6 3 term. 
3. P(s) for Unstable Motion from Perturbation Theory 

Before proceeding further in finding the solutions of equations of motion for unstable 
motion, and their connection with the beta function, (3(s), it will be helpful to examine a 
result for the beta function found using perturbation theory. In reference [1] , the case was 
studied where a particle doing stable motion with the unperturbed tune vq is perturbed 
by a small gradient perturbation which opens up an unstable stopband around vq = q/2, 
where q is some integer. 

In reference [1], the two solutions of the equations of motion inside the stopband are 
found using perturbation theory, and will be denoted here as xi and X2- Using x\ and X2 
one can find the components of M, and in particular 

M\2 = Tj7 {-Xl (s) X 2 (so) + X2 (s) Xl (s )} 

T „ dx 2 dxi 

W = xi— %2—r- 

ds ds 

In the following, the large accelerator approximation is being used, A12 = 1, An = A22 = 0. 
W corresponds to the Wronskian and is a constant of the motion and can be evaluated at 
any value of s. One can find j3 from Myi evaluated at s = sq + L, and Eq. (2-9), 

(3 = (-l) g Mi 2 /sinh^ (3-2) 

The result for xi to lowest order [1], is 

xi = pl (s) exp (g9) cos (q9/2 - (81 + 6 2 ) /2) 



9 = 



\Au\ 2 - (q/2 - z/ ) 2 



A " = hJ dS/3 ° {S) 6XP (~ iq6) 
d9 = ds /vofio, Si = ph (A^) 

5 2 = ph [(q/2 - uq) + ig] 
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The gradient perturbation is AB y = —G(s)x. g is positive. The result given by Eq. (3-3) 
has an error which is first order in the perturbation. The X2 solution is obtained from Eq. 
(3-3) by replacing g by —g and 82 by — 62- fto, vq are the unperturbed beta function and 
tune. 

One may note that x\ can also be written as 



x x = (3* exp [(q + iq/2) 9} [1 + exp [i {q$ - (8 l + 5 2 ))}} 



(3-3b) 



which makes more evident the floquet form of the solution. One sees that /j,, the phase 
change in one turn is given by 

[i = 2%g + iqn 

One can write x\ and X2 as 



xi = /? 2 exp (g8) fa (9) 



x 2 = ffiexp(-g9) h 2 (9) 
hi (9) = cos (qO/2 - (5i + 8 2 ) /2) 
h 2 (9) = cos (qO/2 - (5i - 8 2 ) /2) 



(3-4) 



One then finds 



dx2 
~d9 

X\dX2 

d9 



= Pq exp (-g6) h 2 
= Pohih2 



-g + h' 2 /h2 + -P , /P0 



g + h' 2 /h 2 + ^P'o/Po 



(3-5) 



W = — [-2ghih 2 + hih' 2 - hfa] 

W = — \-2g cos (q9/2 - (61 + 8 2 ) /2) cos (q9/2 - (61 - 8 2 ) /2) 
^0 

- (q/2) cos (q9/2 - + 8 2 ) /2) sin (g(9 2 - (ft - <5 2 ) /2) 

+ q/2 cos (g0/2 - - 8 2 ) /2) sin (g0/2 - (ft + <5 2 ) /2)] 

W = —^-sin8 2 , sin ^2 = (7/I 
2vq 

where the 2gh\h 2 term in W has dropped as being of high order than the remaining term. 
One finds for M 12 from Eq. (3-1) 

M12 = ^ (A) (#) A) (<?o))' [- exp (9 - (9 )] /ii ((9) h 2 (9 ) + exp (9 - 9 )] h x (9 ) h 2 (9)} 

(3-6) 
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Putting 9 = 9 + 2tv and noting that hi(6) = (-l) q hi(9 ), h 2 (0) = h 2 (9 ), and 

2tt5t = n Rl one gets 

M12 = ^ 2 sinh m cos (g(9 /2 - (<*i + 5 2 ) /2) cos - (<?i - <5 2 )) (3-7) 

P can then be found using Eqs. (3-2), (3-5) and (3-6) 

P (9) = p (9) i^^i cos ( q e/2 - (S 1 + 5 2 ) /2) cos (q9/2 - (S 1 - S 2 ) /2) (3-8a) 
Q 9 

Eq. (3-8a) can also be written as 



(9) = p (9) ^o^d [ cos s 2 + cos (q9 - 



9 f 







(3-8b) 



|Az/| 2 - (g/2 - uqY 

Eqs. (3-8) show that as a function of 9, P will change sign and go through zero twice 
in each interval of 2it/q. These two zeros are located at 

qB = Si ± S 2 + 7T (3-9) 

As a function of vq, P becomes infinite near the edge of the stopband \q/2 — uo\ = \Av\ and 
drops to a value of the order of Pq near the center of the stopband, vq = q/2. In one turn, 
A9 = 2/T, P has 2q zeros. One sees that the number of zeros of P in one turn is related to 
the imaginary part of fi, which is qiv. 

This result that connects the number of zeros in P in one turn with the imaginary part 
of the change in the betatron phase in one turn, was found here using perturbation theory, 
but will be shown to be generally valid in the next section. 

4. Solutions of the Equations of Motion and the Beta Function 

For stable motion, the role of the beta function in the solutions of the equations of motion 
is well known. A similar result will be found here for unstable motion. The treatment 
usually given for stable motion, does not carry over to unstable motion because of the 
assumptions usually made about the realness of the betatron and phase functions, and the 
absence of zeros in the beta function. 

Let us write the solutions of the equations of motion as 

x = 6exp (ifj) 

(4-1) 

6 = /32 
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where (3 and b have been defined by Eq. (2-8). Then b has been shown to obey, see Eq. 
(2-18), 

d% T „ 1 
— +Kb + -, = 

ds 2 b 3 (4-2) 

K = -A 21 

The treatment given in this section is for large accelerator case which assumes An = 
A22 = and A21 = 1. Similar results can be found for the general case, x then obeys the 
equations 

d 2 x 

-^ + Kx = (4-3) 

Putting the form of x assumed in Eq. (4-1) into Eq. (4-3), and using Eq. (4-2) for b one 
gets 

dV 2dbdib fdib\ 2 1 

l^ 2 + bd~s~ds~ + \ds~ ) ~¥~ { ' 

Putting / = difj/ds one gets 
The solutions of Eq. (4-5) are 

/ = ± (1/6 2 ) = ±1/0, (4-6) 

Thus 

* = ± fj (4 " 7) 

and the two solutions of the equations of motion are 

-^""Mf) (4 " 8) 

One may note that in deriving Eq. (2-8) no assumption was made about the realness 
of (3 or ip. However, there is a problem with the result for unstable motion, as in the case 
of unstable motion (3(s) will go through zero. To evaluate the integral when (3(s) has zeros, 
Eq. (4-2) will be replaced by 

ib = lim / (4-9) 

^0 J S0 (3-ie 

where e is a positive small quantity. It can be shown that Eq. (4-9) gives (see section 6) 

f ds ds sr^ , 

* = W j^Wmi (4 - 10) 



SO 
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where s n are the locations of the zeros of (3(s) from so to s. P represents the principle 
part of the integral. 

One can also show that /3'(s) = ±2 at the zeros of /3(s). Since (3^ = a 2 — 1, then 
a = ±1 when (3 = 0. Since j3' = —2a, (3' = =f2 when (3 = 0. One can now write Eq. (4-10) 

as 

where N z is the number zeros in (3(s) in so to s. 

One may notice that the imaginary part of tp has on an unusual dependence on s. It 
is constant in between zeros of (3(s) and jumps by tv/2 at each zero of (3(s). One can use 
Eq. (4-11) to find the change in ip over one turn, ip(s + L) — ip(s), and find 

r s+L ds 

ip(s + L)-il>(s) = P -r + iqit (4-12) 

J s P 

where 2q is the number of zeros in (3(s) in one turn, and L is the length of one turn. 
For simplicity, it is being assumed that the period L is one turn. Since (3(s) is a periodic 
function, the number of zeros of (3(s) in one turn has to be even. If one defines the tune 
as the imaginary part of ip(s + L) — ip(s) divided by 2n, then one has 

tune = q/2. (4-13) 

Eq. (4-13) shows the connection between the tune and the number of zeros in the beta 
function in one turn. The real part of i(j(sq + L) — ip(so) gives the exponential growth in one 
turn. If one defines the exponential growth factor, g, to be the real part of tp(s + L) — ip(s) 
divided by 2tv 



P f S+L(1S (A ~\ A\ 

tj= ^L j (4 - 14) 



Another apparent difficulty with the solutions given by Eq. (4-8) is that at the s value 
where (3(s) is zero, both solutions appear to go to zero being proportional to (3?. This is 
not possible as the x motion which is a linear combination of these two solutions would 
then also have to go to zero at this s value. It will now be shown that one of the solutions 
will not go to zero at the zeros of (3(s). 

Let si be a zero of (3(s). Then near s = s\, b = (3i goes to zero like (s — s\) 2 . However, 
it is shown in section 6, that near s = s\ that ip become infinite like {1/ (3'{s\)) log(s — si). 
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Note that j3'(s\) = ±2, and that 6exp(^) goes like (s — s\)^{s — si) ± 2. Depending on the 
sign of /3'(s), bexp(tp) may or may not go to zero at s = s\. If bexp(tp) does go to zero, 



then bexp(-ip) will not go to zero. Thus one of the two solutions will not go to zero at 
s = S\. 

It is interesting to note that the solutions given by Eq. (4-8) can be chosen to be real. 
Let us start at the s value so which is assumed to be in a region where f3(s) is positive 
and let si be the location of the first zero in (3(s) after s = sq. In the region so to s\, 
the solution bexp(ip) is real, as (3 and ip given by Eq. (4-9) are both real. After s = si, 
(3 becomes negative and b = f3^ becomes pure imaginary. However ip jumps at s = s\ by 
i7r/2. Thus the solution bexp(ip) remains real just after s = s\. One can continue in this 
way through the entire lattice with (3 and ip changing suddenly after each zero of (3(s) so 
as to keep the solutions real. This result is consistent with the result found in reference 
[1], that the eigenvalues and eigenfunctions of the one period transfer matrix are real in a 
linear half integer stopband. 

4.1. Eigenvectors of the Transfer Matrix 

The eigenvectors of M will now be found in terms of a and ip. It will also be shown that 
the eigenvalues are given by exp(±Ai/>), so that Aifj = fx, where Aip = ip(s + L) — ip(s). 
Starting from x = (3? exp(±^), one can find the corresponding p x from Eq. (2-10) 




Px 



/TV (-a ± 1) exp (±ip) 



(4-15) 




The two solutions can then be written as 




exp [-ip] 



(4-16) 



These two solutions are the eigenvectors of M as 



Mx\ = exp (Aip) xi, Mx 2 
Aip = ip (s + L) - ip (s) 



exp (—Aip) X2 



(4-17) 
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and the eigenvalues of M are 

Ai = exp (Aip) A 2 = exp (-A^) (4-18) 

Comparing Eq. (4-18) with Eq. (2-7), one sees that fi = Aip = tf;(s + L) —ip(s). 
Since xj exp(±//s/L) is a periodic function, one can write x as 

x = /3* exp [±/ms/L] f (s) 

x = {& exp [±2n (g + iq/L) s/L] f (s) 

where f(s) is periodic with period L. q here is defined by 2q is the number of zeros in f3(s) 
in one turn. One also has 

H = 2n(g + iq/2) (4-20) 

where it is assumed that the period is one turn. 

To summarize, it has been found that if 2q is the number zeros in the beta function in 
one turn, then the eigenvalues of the one period transfer matrix, exp(±//) are given by 

fx = 2ixg + iq% 

P f L ds (4-21a) 



9 2tt y (3 

where P indicates the principle part of the integral, and the solutions of the equations of 
motion are given by 

x = /?2 exp (±-i/>) 

. .. f s ds (4-21b) 
ip = hm / 

J so P- %e 



5. The Emittance Invariant 

The emittance invariant can be found from the Lagrange invariant for symplectic motion. 
If x\ and X2 are two solutions of the equation of motion then [2] 

X2 Sxi = constant (5-1) 

It is assumed that the lattice is periodic, so that the coefficients in the linearized equations 
ofmotion are periodic in s with the period L. Thus if x(s) is a solution then x(s + L) or 
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Using Eq. (2-6) for M one finds 



x sMx = — sin fx \—^x + 2axp x + pp x ) 



(5-3) 



Thus Eq. (5-3) gives the emittance invariant 



e = 7X — 2axp x — Pp x 



e = i (x 2 - (ax + (3p x ) 2 ^j 
(3 1 = l-a 2 



(5-4) 



Eq. (5-4) shows that the curve e = constant is a hyperbola. In the case of stable motion, 
the curve e = constant is an ellipse and e gives the phase space area enclosed by the ellipse. 
For unstable motion, e does not have a simple interpretation in terms of phase space, also 
e can be negative. 

Eq. (5-4) suggests introducing the new symplectic variable 77, p v where 




(5-5) 



Tf 



[5 2 X , p n = f 3 2 (ax + Pp x ) 



G 



1 



The emittance invariant can then be written as 



2 2 
e = r? -p 



(5-6) 
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5.1. Minimum Amplitude 

Eq. (5-4) shows that the particle will move in a hyperbola. Under certain conditions, the 
particle will first move to smaller x or p x before the amplitude of the motion starts to grow 
exponentially. It will be shown below that x and p x can attain the minimum 

Px,min = (-7e) 3 , 7/3 = 1 -a 2 

Eq. (5-7) shows that x will have a mimimum when (3e > 0, and p x will have a mimimum 
when (3e < if \a\ < 1 or /3e > if \a\ > 1. 

The minimum given by Eq. (5-7) can be computed from Eq. (5-4) by computing de/ds 
and putting de/ds = and either dx/ds = or dp x /ds = 0. 

5.2. Asymptotes and Rotation Angle 

If one plots f3p x versus x, one can ask what are the directions of the symptotes of the 
hyperbola. If S\ and 62 are the angles with the x axis for these asymptotes, then they are 
given by 

tan^i = 1 — a 

(5-8) 

tan#2 = — 1 — a 

These results can be found by assuming the asymptotic expansion for (3p x , f3p x = tan 6 x + 
Co + C-ix^ 1 ... and putting this into the equation of the hyperbola, Eq. (5-4). Collecting 
all the x 2 terms and putting the coefficient of x 2 = gives Eq. (5-8). 

If one plots f3p x versus x, then one can ask through what angle this coordinate system 
has to be rotated to make the hyperbola have its normal form. This rotation angle is given 
by 

2a 
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6. Phase Function Results when (3 has Zeros 

In this section, the result for the phase function, ip, given by Eq. (4-10) will be derived. 
Also, the behavior of tp when s is near the zeros of (3(s) will be studied. 
First, let us consider the case where 

</> = lim f (6-1) 

J SQ (3-ie 

e > 0, and one assumes there is only one zero for j3{s) at s = s\ between s = so to s = s. 
Then, one can write 

, „ r s ds r si+s ds . . 

J sa (3 J Sl s p-it 

where 5-^0 but 5 ^> e. P stands for the principle part of the integral. Near s\ one can 
write (3 = (3'(si)(s — si) + ... and find 

ds r i+s ds 



J ai S P J 8l -S P'(si) 



(s — si) — ie 



<™ -2 , -2 » s = s-si, e = e/p (si) 
(si) 7-«5 + (6-3) 
1 ze 

i 



-Z7T 



If there are many zeros between so to s at s = s n one then finds 

Now, it will be shown that near a zero of (3(s), like s = s\, ip becomes infinite like 

ip ~ ±^log(s - si) (6-5) 



The ± corresponds to the sign of (3'(si). We write tp as 

ip = A + B 



A= ds 

J s a 

B = 



< 1 1 



(3 -ie (3' (si) (s - si) - ie J (6-6) 
1 



-so /3' (si) (s - «i) - ie 
where s is assumed to be close to si but s > si. 



16 



The integral of A has no pole near s = s\ and A does not become infinite at s = s\. 
B can be written as 

f3' (si) J (s - si) - ie 

B = {log (s - si - ie) - log (s -si- ie)} (6-7) 

P (si) 

B ~ ±-log(s - si) 

where, in the last result, only the dominant term that becomes infinite at s = si, has been 
kept and the result (3'{si) = ±2 has been used. 

Thus near s = si, ip becomes infinite and the dominant term is given by 

^~±^log( S -si) (6-8) 

where the ± is chosen to correspond to the sign of (3'(s{). 
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